from numpy import *
import scipy as Sci
import scipy.linalg

print('first example')
n=10
# 1.1
A=mat(random.rand(n,n))
B=mat(random.rand(n,n))
# 1.2
Cb=multiply(A,B)
C=A*B
#print(C)
#print(Cb)

# 1.3
b=B[:,0]
x=linalg.solve(A,b)
print(linalg.norm(A*x-b))
